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Abstract. We study entropic functionals/fluctuations of the XY chain with Hamiltonian 

7 -J_ ^(2)^(2) A , \ „(3) 



2 

where initially the left {x < 0)/right {x > 0) part of the chain is in thermal equilibrium at inverse temperature Pi /fir. 
The temperature differential results in a non-trivial energy/entropy flux across the chain. The Evans-Searles (ES) 
entropic functional describes fluctuations of the flux observable with respect to the initial state while the Gallavotti- 
Cohen (GC) functional describes these fluctuations with respect to the steady state (NESS) the chain reaches in the 
large time limit. We also consider the full counting statistics (PCS) of the energy/entropy flux associated to a repeated 
measurement protocol, the variational entropic functional (VAR) that arises as the quantization of the variational char- 
acterization of the classical Evans-Searles functional and a natural class of entropic functional that interpolate be- 
tween PCS and VAR. We compute these functionals in closed form in terms of the scattering data of the Jacobi matrix 

iUx-1 canonically associated to the XY chain. We show that all these functionals are 
identical if and only if h is reflectionless (we call this phenomenon entropic identity). If h is not reflectionless, then the 
ES and GC functionals remain equal but differ from the PCS, VAR and interpolating functionals. Purthermore, in the 
non-reflectionless case, the ES/GC functional does not vanish at a = 1 (i.e., the Kawasaki identity fails) and does not 
have the celebrated a 1 — a symmetry. The PCS, VAR and interpolating functionals always have this symmetry. 
In the Schrodinger case, where Jx = J for all x, the entropic identity leads to some unexpected open problems in the 
spectral theory of one-dimensional discrete Schrodinger operators. 



1 Introduction 



The XY chain is an exactly solvable model which is often used to illustrate existing theories and to test 
emerging theories in quantum statistical mechanics. The literature on XY chains is enormous and we 
mention here the mathematically rigorous works [Ar3, Ar4, ABl, AB2, AH, AP, BM, HR, Ma, Mc, OM]. 
It is therefore not surprising that a model of an open quantum system based on the XY chain was also a 
testing ground for the recent developments in non-equilibrium quantum statistical mechanics initiated in 
the works [JPl, JP2, Ru2]. More precisely, the first proofs of the existence of a non-equilibrium steady 
state (NESS) and of strict positivity of entropy production in an open quantum system were given in the 
papers [AH, AP] in the context of XY chains. The purpose of this paper is similar and our goal is to test the 
emerging theory of entropic fluctuations in non-equilibrium quantum statistical mechanics developed in 
[JOPP] (see also [DDM, Ku, Ro, TM] for related works) on the example of XY chains. For discussion and 
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references to the classical theory of entropic fluctuations (which, starting with the seminal works [CG, ES], 
has played a dominant role in the recent theoretical, numerical and experimental advances in classical 
non-equiUbrium statistical mechanics) we refer the reader to the reviews [JPR, RM]. 

The paper is organized as follows. In Section 1.1 we describe the XY chain confined to a finite interval in 
Z, introduce the finite volume and finite time entropic functionals and describe their basic properties. The 
thermodynamic limit of the confined XY chain and its entropic functionals is discussed in Section 1.2. The 
elements of spectral and scattering theory needed to state and prove our results are reviewed in Section 1 .3. 
In Section 1.4 we review the results regarding the existence of the NESS and introduce the finite time 
Gallavotti-Cohen functional. Our main results are stated in Section 1.5. In Section 1.6 we briefly discuss 
these results and in Section 1.7 we comment on reflectionless Jacobi matrices. The proofs are given in 
Section 2. Some proofs are only sketched and the reader can find details and additional information in the 
forthcoming Master's thesis [La]. 

Acknowledgment. The research of V.J. and B.L. was partly supported by NSERC. The research of C.-A.R 
was partly supported by ANR (grant 09-BLAN-0098). We wish to thank C. Remling and B. Simon for 
useful discussions. 

1.1 The confined XY chain 

To each a; G Z we associate the Hilbert space Hx = and the corresponding matrix algebra Ox = 
M2(C). The Pauli matrices 

1 

-1 J ' 

together with the identity matrix l^^ form a basis of Ox that satisfies the relations 

Let A = [TV, M] be a finite interval in Z. The Hilbert space and the algebra of observables of the XY chain 
confined to A are 

nA = <^nx, Oa = (^Ox. 

xeA xeA 

The spectrum of the observable A G Oa is denoted sp(A). If A is self-adjoint, 1e{A) denotes the spectral 
projection associated to E E sp{A). 

We shall identify A^ G Ox with the element {^y£A\{x}'^y) ^ of Oa- In a similar way we identify Oa' 
with the appropriate subalgebra of Oa for A' C A. With these notational conventions the Hamiltonian of 
the XY chain confined to A is 

x£[N,M[ xe[N.AI] 

Here, Jx is the nearest neighbor coupling constant and A^; is the strength of an external magnetic field in 
direction (3) at the site x. Throughout the paper we shall assume that {Jx}xez and {Xx}xei, arc bounded 
sequences of real numbers and that Jx for all x G Z. 

Consider the XY chain confined to the interval A = [—M, M]. Its Hamiltonian can be written as 

HA=Hi+Hr + V, 

where Hi is the Hamiltonian of the XY chain confined to [—M, 0], Hr is the Hamiltonian of the XY chain 
confined to [1, and 

V T (2)\ 



1 

1 



r(2) = 



-i 

1 



.(3) = 
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is the coupling energy. Until the end of this section A = [—A/, M] is fixed and we shall omit the respective 
subscript, i.e., we write T-L ~ Ha, O = 0\, H = Ha, etc' 

The states of XY chain are described by density matrices on H. The expectation value of the observable A 
w.r.t. the state p is p{A) = tj:{pA) . We set 

t\A) = e'*" Ae"''" . 

In the Heisenberg/Schrodinger picture the observables/states evolve in time as 

At = T\A), Pt^T-\p). 

Obviously, pt{A) ~ p{At). The relative entropy of the state p w.rt. the state v is defined by 

S{p\v) =tr(/3(logi/-logp)). 

We recall that S{p\v) < and that S{p\i') = iff p = i^.- For additional information about the relative 
entropy we refer the reader to [JOPP, OP]. 

Following the setup of [AH], we shall consider the evolution of the XY chain with an initial state given by 



(1.1) 



where > 0. Hence, the left part of the chain is initially in thermal equilibrium at inverse temperature 
Pi while the right part of the chain is in thermal equilibrium at inverse temperature f3r. The XY chain is 
time-reversal invariant: there exists an anti-unitary involution 9 : H ^ H, described in Section 2.1, such 
that, setting Q{A) = 9A6^^, one has 

e(i/,/,) = Hi/,, Q[V) = V, = u. 

The observables describing the heat fluxes out of the left/right chain are 

d 



and an easy computation gives 



* - /- 7- ^^3), (1) (2) (2) (1)^ 1 . , /^(2)^(1) ^(l)^(2)^ 

The entropy production observable of the XY chain is 

Note that $; and a change sign under the time reversal: 

e($,/,) = e(a) = -a. 



'One can consider a more general setup where Hi is the Hamiltonian of the XY chain confined to [— M, — — l]/[N + l, M], 
He is the Hamiltonian of the XY chain confined to [— A'^, A'^], and H = Hi + He + Hr + V, where 

1/ 1 r / (1) (1) I (2) (2) \ 1 , / (1) (1) I (2) (2) \ 

V = -J-N-1 [^Ik-l'^-N + <^-iv-l<^-JV ) + [(^N^N+l + < ) ■ 

In this scenario the left and the right parts of the chain are connected via the central part which is also an XY chain. In the thermody- 
namic limit one takes AI oo while keeping A'^ fixed. The central part remains finite and represents a small system (quantum dot) 
connecting two infinitely extended chains (thermal reservoirs). All our results and proofs extend to this setting [L,;]. For notational 
simplicity, we will discuss here only the simplest case of two directly coupled chains. 
" S{p\u) = — oo unless Ker u C Kcr p. 
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The observable describing the mean entropy production rate over the time interval [0, t] is 

1 /■* 
S* = - / asds. 
* Jo 

The basic properties of this observable are summarized in: 

Proposition 1.1 (1) logwt = logw + ir^*(E*). 

(2) SiutH = 

(3) r*(I]-*) = S*. 

(4) E* ~ — r*(0(E*)). /« particular, sp(E*) is symmetric w.r.t. zero and dim 10(S*) = dim l_0(S*)/or 
fl// e sp(I]*). 

The proof of Proposition 1.1 is elementary and can be found in [JOPF]. Part (1) allows for the interpretation 
of the entropy production observable cr as the quantum phase space contraction rate. Part (2) implies that 
for all i > the average entropy production rate over the interval [0, t] is non-negative, i.e., 

04 >0, (1.2) 

0esp(S') 

where 

is the probability that a measurement of E* in the state uj will yield the value (p. In particular, on average 
heat flows from the hotter to the colder part of the chain, in accordance with the (finite time) second law of 
thermodynamics. By Property (4), (1.2) is equivalent to 

J2 0(p;-p%)>o. 

0esp(s') 
4>>o 

The finite time fluctuation relation is deeper The direct quantization of the Evans-Searles fluctuation 
relation in classical non-equilibrium statistical mechanics is the identity 

P%=e-*^4, (1.3) 

which should hold for all (f) and t > 0. An equivalent formulation of ( 1.3) is that the functional 

ESt{a) = lege, (c-"*^') = log e-"*V^, 

0Gsp(E*) 

has the symmetry 

ESt(a) =ESt(l-a), (1.4) 

that holds for all a e R and t > 0. The functional ESi(a) is the direct quantization of the finite time 
Evans-Searles functional in non-equilibrium classical statistical mechanics [''^R]. It is however easy to 
show that ESt(l) > = ESt(O) except at possibly countably many <'s and so the relations (1.3) and (1.4) 
cannot hold for alH > (see Exercise 3.3. in [JOPP]). This point is further discussed in [La]. 

The first class of quantum entropic functionals that satisfy the Evans-Searles fluctuation relation for all 
times was proposed independently by Kurchan [Ku] and Tasaki-Matsui [TM]. They involve the funda- 
mental concept of Full Counting Statistics associated to the repeated measurement protocol of the en- 
ergy/entropy flow [LL]. Consider the observable 

5 = - logCJ = PiHi + prHr + Z, 
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where Z =■ logtr(e ^' ' Pr-Hr-y Clearly, St = t (S) = — log cj_t, and Properties (1) and (3) imply 

S* = i(5,-5). 

The probability that a measurement of S at time t = (when the system is in the state a;) yields s € sp(5) 
is a;(ls(iS)). After the measurement, the system is in the reduced state 

UjlsjS) 

which evolves in time as 

A second measurement of 5 at a later time t > yields s' E sp(5) with probability 

and the joint probability of these two measurements is 

tr (c-"^wl,(5)c"^l,,(5)) . 

The mean rate of entropy change between the two measurements is (p = {s' — s)/t and its probability 
distribution is 

Pt(</))= tr(e-'*^c^l,(5)e'*^l,-(5)). 

s' —S—t(p 

The discrete probability measure Pj is the Full Counting Statistics (PCS) for the operationally defined 
entropy change over the time interval [0, t] as specified by the above measurement protocol. Let 

FCSt(a) -log^e-"*^Pt(,^). 

One easily verifies the identity 

FCSt(a) =logtr(c.i-"w"), (1.5) 
and time-reversal invariance implies the fluctuation relation 

FCSt(a) =FCSt(l-a), (1.6) 

which holds for all a and t. The last relation is equivalent to 

Pt(-0) =e-*'^Pt ((/)). (1.7) 

The identities (1.6) and (1.7) provide a physically and mathematically appealing extension of the Evans- 
Searles fluctuation relation to the quantum domain. 

Generalizations of the fluctuation relation (1.6) have been recently proposed in [JOPP]. Note that 

FCSt(a) = logtr (e^-^e^-^'e^ 
Forp > and a e M we define the functionals ((a) by 



ep,t{.a) 



2°5j l^s 

CP e p e p 



if < p < oo, 



log tr 

logtr (e^i-^^-^+aSt-) ifp^oo. 
Their basic properties are summarized in; 
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Proposition 1.2 (1) ep,t(0) = ep,t(l) = 0. 

(2) The functions W 9 a i— > ep^t(a) are real-analytic and convex. 

(3) ep^t{a) = ep^-t{a). 

(4) ep.tia) = ep,t(l - a). 

(5) The function ]0, cx)] 9 p i— > ep.t(a) continuous and decreasing. 

(6) Cp ((0) = ESj(O) = — j(l) = — ta;(E*). /n particular, these derivatives do not depend on p. 

(7) e2^t(a) =FCSt(a) W 

e^'_,(0) = ES;'(0)= / / -cj(a,))(a„-cjK)))dsdw. 
Jo Jo 

(8) Boo, t(ct) = max {S{p\ijj) — atp(Y:*)), where the maximum is taken over all states p. 

Remark. The variational characterization of etx3,t(a) in (8) is the quantization of the variational char- 
acterization of the finite time Evans-Searles functional in classical non-equilibrium statistical mechanics. 
Regarding the physical interpretation of ep_t(a) forp < oo, note that 



tr (^c^^c^^'e^^^ /tr(w2/p) ^ J e-"*'^dPp,t (</>), 
where Pp t is the PCS of the XY chain with scaled temperatures 2/3i/r/p- Restoring the scaling as 



tr 



^St 

e p e p c p 



p/2 



distorts the connection with the PCS but links ep^t{a) with quantized Ruelle transfer operators which 
provide yet another way of quantizing the classical Evans-Searles functional (see [JOPP, JOP, JP3]). 

Por the proof of Proposition 1.2 we refer the reader to [JOPP]. The definitions and structural relations 
described in this section have a simple general algebraic origin and are easily extended to any time-reversal 
invariant finite dimensional quantum system [JOPP]. A similar remark applies to the thermodynamic limit 
results described in the next section which hold for a much wider class of models than XY chains and in 
particular for quite general lattice quantum spin systems (for example, for the open spin systems discussed 
in [Ru2]). However, our study of the large time limit t — > od critically depends on the specific properties 
of the XY model. 



1.2 Thermodynamic limit 

In this section we describe the thermodynamic limit A = [—M, M] /• Z (for discussion of the thermody- 
namic limit of general spin systems we refer the reader to [BRl , BR2, 1, Ru 1 , S]). We use the subscript M 
to denote the dependence of various objects on the size of A and write, for example, Om, Hm, ep,t,M, 
etc. The C* -algebra O of observables of the extended XY chain is the uniform closure of the algebra of 
local observables 

M 

where we identify Omi with a subalgebra of Om^ for Mi < A/2- Por any A G 0\oc the limit 

t\A) = lim c''"" Ac~'*"" , 

exists in norm and t* uniquely extends to a strongly continuous group of *-automorphisms of O that 
describes the dynamics of the extended XY chain. The physical states of the extended XY chain are 
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described by positive normalized linear functionals on O. The expectation of an observable A in the state 
p is p{A). In the Heisenberg/Schrodinger picture the observables/states evolve in time as 

At=T*{A), pt^porK 

S{p\i') denotes Araki's relative entropy of the state p w.rt. the state v [Arl, Ar2] with the notational 
convention of [BR2]. Let unj be the reference state (1. 1) on Om- Then for all A g Oioc the limit 

uj{A) = lim ujm{A), 

M^oc 

exists and lu uniquely extends to a state on O that describes the initial state of the extended XY chain. 
The C* dynamical system (O, r*,a;) describes the extended XY chain. This quantum dynamical system 
is time-reversal invariant, i.e., there exists an anti-linear involutive ^-automorphism 8:0^0 such that 
e o r* = r^* o e for all t and w(e(A)) = lo{A*) for all AeO. 

The observables and a are obviously in Oioc- Let 

1 /■* 
E* = - / asds, 

t Jo 

ESt(a) = logw(e-"*^'). 

We then have 
Proposition 1.3 (1) 

lim E*M = S*, 

A/— i-oo 

in norm and the convergence is uniform for t's in compact sets. 

(2) S{LOt\Lo) = -tuj{T.^). 

(3) r*(I]-*) ^ S*. 

(4) = — T*(0(I]*)). In particular, sp(I]*) is symmetric with respect to the origin. 
(5) 

lim ESf.A/(a) ESt(a), 
and the convergence is uniform for t's and a's in compact sets. 

Regarding the entropic functionals, we have: 

Proposition 1.4 (1) For alia £R andp e]0, oo] the limits 

ep,tia) = lim ep^t.Aiia), 

M—^oo 

exist and are finite. 

(2) ep,t(0) = ep,t(l) =0. 

(3) The function K 9 a k->- ep^t{ct) ore real-analytic and convex. 

(4) Cp^tia) = ep-t{a). 

(5) ep_i(a) = ep,t(l - a). 

(6) The function ]0, cx)] 9 p i— > ep^i(Q;) w continuous and decreasing. 

(7) e;,(0) = ES;(0) = -e;,,(l) = -to.(E*). 
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(8) 

e'lt{0) = ES"(0) = / / cj {{as - w(crs))(cr„ - w(cr„))) dsdu. 
Jo Jo 

(9) As M — >■ oo, the sequence Pj^a/ converges weakly towards a Borel probability measure Pt on K. and 

e2,t(a) = log / e-*"^dPt ((/)). 

All the moments of¥t^M converge to corresponding moments o/P^. The measure P^ is the full counting 
statistics of the extended XY chain. 

(10) eao,t{oi) = sup {S{p\uj) — atp{Yi^)), where the supremum is taken over all states p. 

The results of this section are easy to prove (see Section 2.3 and [La]). Also, using standard techniques 
[BRl, BR2, 1, Rul, S] these results are easily extended to other models in quantum statistical mechanics 
(general lattice quantum spin system, interacting fermionic models, Pauli-Fierz models, etc.). We empha- 
size, however, that non-trivial thermodynamic al behavior only emerges in the large time limit t ^ oo and 
that controlling the large time limit of entropic functionals in physically interesting models is typically a 
very hard analytic problem (see [Ro, JOPP] for existing results and discussion of this point). The advantage 
of models like the XY chain is that the large time limit is easy to control via trace class scattering theory 
and very detailed information about limiting entropic functionals is available. In the next section we review 
some basic results of spectral and scattering theory that we will need in this paper 

1.3 Operator theory preliminaries 

Let A be a self-adjoint operator on a Hilbert space S^. We denote by spa(,(A) the absolutely continuous 
spectrum of A and by i3ac(^) the corresponding spectral subspace. The projection on this subspace is 
denoted by lac(^). For any ipi,ip2 S i5 the boundary values 

(^1, {A-E- iO)- V2) = lim(V;i, {A-E- ie)" V2), 

exist and are finite for Lebesgue a.e. ii^ € M. In what follows, whenever we write (-01, {A — E — i0)~^ip2), 
we will always assume that the limit exists and is finite. If ipi = il;2 = ip then Im (?/', (A— i?— iO)^^i/') ^ 0- 
If i/^ is the spectral measure of A for ip, then its absolutely continuous component with respect to Lebesgue 
measure is 

diy^ME) = -Im {'^, {A-E- iO)-^^)dE. (1.8) 
For the proofs and references regarding the above results we refer the reader to [J]. 

The basic tool in virtually any study of XY chains is the Jordan- Wigner transformation [JW] (see also 
[LSM, Ai'3, AP, JOPP]). This transformation associates to our XY chain the Jacobi matrix 

{hu){x) = Jxu{x + 1) + \xu{x) + Jx-iu{x — 1), 

(see Section 2.1). As an operator on the Hilbert space [) = i^{'L), h is bounded and self-adjoint. Set 

i)i^e\]-cx),o]), i)r^e^{[i,^[), 

and denote hi/^ the restrictions of h to f);/^ with a Dirichlet boundary condition. Note that h = ho + v 
where ho = hi (B hr and 

V = Jo{\Sl){Sr\ + \dr){Sl\), (1.9) 

61 /r being the Kronecker delta function at site x — 0/1. We denote by vi/r the spectral measure of hi/^ for 
,5,/,. By (1.8), 

d^^//,,ac(S) = -Fi^,.{E)dE, 
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where 

Fi^r{E) = ImGifriE), Gi/r{E) - {5i/,., {hi,, -E- i^rHi,,). 

By the spectral theorem, one can identify f)ac(/io) with L^(R, di^;_ac) ffi i^(M, dvr^ac) and /lo |" f)ac(^o) 
with the operator of multipHcation by the variable £' e M. The set S; /^ ac = G I > 0} is the 

essential support of the absolutely continuous spectrum of hij^. We set 

^ — ^/,ac ri Xj^ ac- 

Finally, we recall a few basic facts that follow from trace class scattering theory [RS, Y]. The wave 
operators 

zz;± =s-lime'*''c-"''"lac(/io), 

i— ^±00 

exist. The scattering matrix 

s = 'wXw- (1.10) 

is unitary on \)s.c{ho) and acts as an operator of multiplication by a unitary 2x2 matrix function (called 
the on-shell scattering matrix) 

'^""^ " [ Srl{E) Srr{E) 

where 



Sab{E) = 5ab + 2\4{5a, {h - E - 5^) ^ Fa{E)Fb{E) . 

In our current setting these results can be easily proven directly (see [JKP, La]). Note that s{E) is diagonal 
for Lebesgue a.e. £^ e M \ It follows from the formula 

{5u {h-E~ iO)-%) = {5r, {h-E~ iOr'di) = - '^"^'(^)^^(^) 



I - J§GiiE)GriE) 
that s{E) is symmetric and not diagonal for Lebesgue a.e. E £ £. 

The Jacobi matrix h is called reflectionless iff s{E) is off-diagonal for E E £. In other words, h is 
reflectionless if the transmission probability satisfies \sir{E)\'^ = 1 for E E £. For other equivalent 
definitions of reflectionless we refer the reader to Chapter 8 in [Te] (see also [La] for a discussion). For 
additional information and references about reflectionless Jacobi matrices we refer the reader to the recent 
work [Re]. 

1.4 The large time limit: NESS 

The basic result concerning the existence of a non-equilibrium steady state (NESS) of the XY chain is: 

Theorem 1.5 Suppose that h has purely absolutely continuous spectrum. Then for all A E O the limit 

(A)+= lim w(t*(A)), 

exists. The state Lu^{-) = {■)-^- is called the NESS of the quantum dynamical system {O^t* ,lu). The 
steady state heat fluxes are 

1 /■ , , ,,n smh(ABE/2) 
- = y,^'--'-(^)' cosh(/3.V2)cosh(Ai./2) ^^' ^'-''^ 

and the steady state entropy production is 
where /\p = fir - Pi. 
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With only notational changes the proof in [AP] extend to the proof of Theorem 1.5\ If in addition to its 
absolutely continuous spectrum h has non-empty pure point spectrum then one can show that the limit 



{A)+ = lim - / uj{T'{A))ds, 



T-^'^ T Jo 

exists for all A E O and that the formula (1.11) remains valid (see [AJPP2]). In presence of singular 
continuous spectrum the existence of a NESS is not known. 

If \£\, the Lebesgue measure of £, is zero, then obviously there is no energy transfer between the left and 
the right part of the chain and {a)+ = ($i/r)+ = for all (3i/r. If \£\ > 0, then {a)+ > iff f3i ^ Pr. 
Moreover, ($;)+ > if /?; < /3r and (^>;)+ < if /3i > (3r- The formula (1.11) is of course a special case 
of the celebrated Landauer-Biittiker formula that expresses currents in terms of the scattering data."* 

Fluctuations of the mean entropy production rate S* with respect to the NESS u;+ are controlled by the 
functional 

GCt(a)=log^+(c-"*^'), 

which is the direct quantization of the Gallavotti-Cohen functional in non-equilibrium classical statistical 
mechanics [JPR]. 

1.5 The large time limit: Entropic functionals 

In this section we state our main results. It is convenient to introduce the following matrix notation. We let 



-piE 

-PrE 



set 

Ka{E) = Qko(E)/2^a(s'(E)ko(E)s{E)-ka(E))^ko(E)/2 ^ ^.^^2) 

and note that Kq{E) = c^o^^) . We further set 

p/2 



(1.13) 



forp g]0, oo[, and 

^a,oo(£^) = lim K^JE) = Q(^-°')k«(E)+o.^{E)ko(E)s'{E) _ 
' p— >oo 

Note that Ko.p{E) = Ko{E). Moreover, it is a simple matter to check that, for E,a and p g]0, oo], 

Ka{E) = Ka{E) ^ Ka^p{E) = Ko{E) ^ [s{E), ko{E)] = 4^ A = or s{E) is diagonal. 
Recall that the last condition holds for Lebesgue a.e. E gM.\£, but fails for Lebesgue a.e. E E £. 

Theorem 1.6 Suppose that h has purely absolutely continuous spectrutn. Then the following holds: 

(1) For a e R and p g]0, oo], 

.^ 1 /, f dct{l + K^,piE)) \ dE 



lim -ESt(a) = lim -GCt(a) = e+(a) = / log - — 7 ' -— . 

These functionals are identical to zero iff\£\ = Q or Pi ~ /3r. In what follows we assume that |£| > 
and Pi ^ Pr- 



^Alternatively, Theorem 1.5 can established by applying the results of [JKI^] to the Jordan-Wigner transformed XY chain. 
We refer the reader to [/\JPP2, N, IKP] for mathematically rigorous results regarding the Landauer-Biittiker formula and for 



references to the vast physical literature on the subject. 
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(2) ThefunctionM. 9 a i— > ep_+(a) is real -analytic and strictly convex. Moreover, ep^+(0) ~ 0, _|_(0) = 
— (cr)_l-, and 

ep,+(a) = ep,+ (l - a). 

(3) The function M 9 a i— > 6+ (a) real-analytic and strictly convex. Moreover, it satisfies e+{Q) = 0, 
e+(0) = — (ct)+, and 

e'KO) = e'{+(0) = ^lirn^ ^ j\ ^'^"^ ' + " + '^^^ ^^'^^^ 

(4) e+(l) > unless h is reflectionless. Ifh is reflectionless then 

/cosh((;3/(l - a) + Pra)E/2) cosh((/3^(l - a) + Pia)E/2)\ AE 
\ cosh(/3,£;/2)cosh(^,.i;/2) J 2^' 

flnc/e+(a) = e-|_(l — a). 

(5) The fiinction ]0,oo] 3 p ep -|_(a) /i continuous and decreasing. It is strictly decreasing for a ^ 
{0, 1} unless h is reflectionless. Ifh is reflectionless, then ep_+(a) does not depend on p and is equal 
to e+(a). 

Remark 1. Typically, the correlation function C{s) = {{as — {a) + ){a — (<t) + ))^ is not absolutely in- 
tegrable and the Cesaro sum in (1.15) cannot be replaced by the integral of C{s) over R without further 
assumptions. 

Remark 2. As already noticed, time-reversal invariance implies that the on-shell scattering matrix is 
symmetric. It follows that one can replace s{E) by s*{E) and s*{E) by s{E) in (1.12)-(1.14) when 
inserted in the formulas of Part (1 ). 

Remark 3. The vanishing of entropic functionals at a = 1 is called Kawasaki's identity, see [CWWSE]. 
Thus, the Kawasaki identity holds for all e^'s and fails for e+ unless h is reflectionless. 

Theorem 1.6 has the following implications. To avoid discussing trivialities until the end of this section we 
assume that |£| > and /3; ^ (3,.. Recall that 

e2,t(a) = FCSt(a) - / c'^^'^dPt ((/.), 

where Pt is the PCS measure of the extended XY chain. If Pes,* and Pact are respectively the spectral 
measures of E* for lu and w+ we also have 

ESt(a) = log / b-'^'UPbsA^), GCtia) = log / e-"*'^dPGc,t(0). 
The respective large deviation rate functions are given by 

Ifcs+{(}) = - inf {ad + e2.+ {a)), 
I+{9) = -mUa9 + e+{a)). 

The functions /+ [9) and /fcs+ {&) are non-negative, real-analytic, strictly convex and vanish at the single 
point 9 = (cr) + . These two rate functions are different unless h is reflectionless. The symmetry 62,+ (a) = 
62, +(1 — a) imphes 

^FCS+(^) = ^FCS+( — + ^• 

The rate function I+{9) satisfies this relation only if h is reflectionless. Theorem 1.6 implies: 
Corollary 1.7 Suppose that h has purely absolutely continuous spectrum. 
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(1) The Large Deviation Principle holds: for any open set O C M, 

lim ilogPEs,t(0) = lini ilogPGC,t(0) - - inf 

t-i-oo t t->oo t 9eo 

lim ilogPFcs.t(O) = - inf /fcs+(^). 

(2) The Central Limit Theorem holds: for any Borel set S C R, let Bt ^ y/t{(p - (<t) + ) G B}. Then 

lim PEsABt) = lim PGC,t(-Bt) = lim FpcsABt) = J— [ e""^' /^""^ d^, 

t^oo t->-oo t-^oo ^2txBj^ Jb 

where the variance is = e" (0). 

1.6 Remarks 

Due to non-commutativity it is natural that non-equilibrium quantum statistical mechanics has a richer 
mathematical structure than its classical counterpart. This is partly reflected in the emergence of novel en- 
tropic functionals. The direct quantizations of the Evans-Searles and Gallavotti-Cohen entropic functionals 
typically will not have the symmetries which, starting with the seminal works [CG, ES], have played a 
central role in recent developments in non-equilibrium classical statistical mechanics. The mathematical 
theory of entropic fluctuations in quantum statistical mechanics is an emerging research direction and our 
testing of the existing structural results on the specific example of the XY chain has lead to some surprising 
(at least to us) results. 

If in a given model all the functionals ep,+ {a) are equal and coincide with e+(a) we shall say that in this 
model the entropic identity holds. For the XY chain we have shown that the entropic identity holds iff 
the underlying Jacobi matrix is reflectionless. Note that the large time fluctuations of Pss.t and ^GC,t are 
identical for all XY chains (i.e., for any h) and this is not surprising — the same phenomenon generally 
holds in classical/quantum non-equilibrium statistical mechanics (this is the principle of regular entropic 
fluctuations of [JPR, JOP]). On the other hand, the full counting statistics of repeated measurement of 
entropy flux is obviously of quantum origin with no classical counterpart and in general one certainly does 
not expect that the large time fluctuations of Ppcs,* are equal to those of Pes/GC,*- In the same vein, one 
certainly expects that eoo,+ (a) < 62,+ (a) for a 7^ 0, 1, and more generally, that all functionals ep^+{a) 
are different. 

Needless to say, there are very few models for which the existence of ep +(0;) and e+{a) can be rigorously 
proven and the phenomenon of entropic identity, to the best of our knowledge, has not been previously 
discussed in mathematically rigorous literature on the subject. The XY chain is instructive since the func- 
tionals ep^+{a) and e+(a) can be computed in closed form and the entropic identity of the model can be 
identified with the reflectionless property of the underlying Jacobi matrix. The question whether such a 
physically and mathematically natural result can be extended beyond exactly solvable models like the XY 
chain or the Electronic Black Box Model^ remains to be studied in the future. 

1.7 Examples 

Schrodinger case. If ~ J for all x then the Jacobi matrix h is the discrete Schrodinger operator 

{hu)(x) = J{u{x + 1) + u{x — 1)) + \xu{x). 

In the literature, the most commonly studied XY chains correspond to this case. The assumption that h has 
purely absolutely continuous spectrum plays a critical role in the formulation and the proof of Theorem 
1.6. The striking fact is that the only known examples of discrete d — 1 Schrodinger operators with 

^See [AJJPl, AJPP2, JKP, JOPP] for discussion of these models and references. 
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purely absolutely continuous spectrum are reflectionless. If the potential is constant or periodic then h 
has purely absolutely continuous spectrum (and is reflectionless). The only other known examples involve 
quasi-periodic potentials. For example, if A^; — cos(27rQ:a;+0), then for |J| > 1 /2 the operator /t, has purely 
absolutely continuous spectrum and is reflectionless (see [A\ ]). In this context it is an important open 
problem whether there exists a one-dimensional discrete Schrodinger operator h which is reflectionless 
and has purely absolutely continuous spectrum and therefore an XY chain with constant for which the 
entropic identity fails 

Jacobi case. If the are allowed to vary, then it is easy to produce examples of h which are not reflec- 
tionless and have purely absolutely continuous spectrum (take A^; = 0, J^; = 1 for .t > and = 1/2 
for X < 0). In fact, the vast majority of Jacobi matrices with purely absolutely continuous spectrum are not 
reflectionless. To illustrate this point, let q : [—2, 2] C be such that < 1 and 

r 'J^aE < oo. 

J-2 V2^2 

Then there exists a Jacobi matrix h such that the spectrum of h is purely absolutely continuous and equal to 
[—2, 2], and such that sir{E) = q{E) for Lebesgue a.e. E G [—2, 2] [VY]. The only reflectionless Jacobi 
matrix in this class is the Schrodinger operator with J = ±1 and A^; = (see [Te]). A similar result holds 
if the interval [—2,2] is replaced by a homogeneous set (say, a finite union of intervals or even a Cantor set 
of positive measure) [VY]. 



2 Proofs 

2.1 The Jordan- Wigner transformation 

Consider the XY chain confined to A = [— M, M] . In this subsection M is fixed and we drop the respective 
subscript. Let A( = 0], A^ = Set fi = ^^(A), = £'^{Ai/r), and let h, hi/r be the 

restrictions of the Jacobi matrix to A, A;/^. Whenever the meaning is clear we extend hi/j. to fi by ft,/ © 
0/0 ® hr- Let ho = hi + ft,, and 

ft = fto + 

where v is given by (1.9). 

We shall assume that the reader is familiar with the formalism of the fermionic second quantization (see 
[BR2, BSZ, AJJPl] for general results and [JOPP] for a pedagogical introduction to this topic). T denotes 
the antisymmetric Fock space over i) and a* {f)/a{f) the creation/annihilation operator associated to / € f) . 
We write a*{Sx) ~ a* etc. Set S-m — 1 and 

5,= n (2a;a,-l), 

ye[~M,x[ 

foTx e] ~ M,M]. 

Proposition 2.1 There exists a unitary ?7jw : H — > J^, called the Jordan-Wigner transformation, such 
that 

UjwCTx^Uj'^ = Sx{ax + a*), Ujvfa^^^UJ'^ = iS^^ia^ - a*), UjwcTx^Uj'^ ^ 2a*ax - 1. 

The Jordan-Wigner transformation goes back to [J W] and the proof of the above proposition is well-known 
(a pedagogical exposition can be found in [JOPP]). An immediate consequence are the identities 

UjwHifrUj-^ - dr(ft,/,), UjwVUj-^ = dr(i;), UjwHUj-^ = dr(ft), 

*It is believed that such an example exists. We are grateful to C. Remling and B. Simon for discussions on reflectionless Jacobi 
matrices. 
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Ujw^iUj^ = -iJo J-i (a*a_i - a*_^ai) - UqXo (a^ao - a^ai) , 
J7jw$rC^j^ = -i-'^o-^i (ata2 - a^ai) - Uo\i {a^ai - alao) . 

It is also easy to see that the confined XY chain is time-reversal invariant with 9 = C/j^r(j)J7jw, where 
ji/; = tp is the standard complex conjugation on £^(A). 

In what follows we will work only in the fermionic representation of the confined XY chain. With a slight 
abuse of notation we write Hi/^ = dT{hi/j,), V = dr(ii), H = Hi + Hr + V = dT{h), etc. 

2.2 Basic formulas 

In this section A = [—A/, M] is again fixed and we drop the respective subscript. We shall make repeated 
use of the identity 

tr(r(A)) =^ dct(l + A), 
which holds for any linear map A : t) ^ t). Set 

k = ~l3ihi - PrK, h = e'*''fce-'*''. 

The initial state of the system is described by the density matrix 

r(e'=) 



tr(e-A«'"'3'-^0 det(l + efc)' 



Since 



we have 



V * / dct(l + c'=) 



ep,t(a) = loj 



dct(] 

det (l + (e(i-")fe/Pe2"*^--*/f e(i-«)'=/P)''/') 



det(l+e'=) 
for < p < oo. Similarly, 

det (1 + 6(1-")'=+"'=-') 
^-'^")=^°^' det(l+e^) ' 

and 

det('l + e'=/2e«('=*-'=)e'=/2) 
^ det(l + e^) 

Clearly, the functions a ^ ep^t {a), a ESt(a) are real analytic. For u S R and < p < oo set 

/Cp,t(a,w) = ie-(i-")'=*"/P ^1+ (^e(i-")'='"/Pe2"'=-'<i-")/Pe(i-")'='"/p) ''^^^ e^i-")'^'"/^ + h.c, 

(2.16) 

where h.c. stands for the hermitian conjugate of the first term. We also set 

/Coo,t(a,u) = (l + c-(i-")'=--"'=-*(i-"))"\ (2.17) 



and 

/CES,t(a,M) = - (l + e-"('='(i-")-'=-'")e-'=-'") . (2.18) 
The following lemma will play the key role in the sequel. 
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Lemma 2.2 (1) Forpe]0,oo], 

ep,t(a) = i/ d7 / dutr(/Cp,t(7,u)i[A:,/i]). (2.19) 
Jvi Jo 

(2) 

pa. pi 

ESt{a)^t d7 / dutr(/CEs,t(7,")i[fc,/^])- (2-20) 
Jo Jo 

To prove the lemma we need the following preliminary result (for the proof see Lemma 2.2. in [HP]): 

Lemma 2.3 Let F be a dijferentiable function of the real variable a taking values in Hermitian strictly 
positive matrices on f). Then for any p €]0, oo[, 

^iv\og{l + F{aY)=piT ( {1 + F{a)-P)-^F{a)-^^^^ 
da \ da 

Proof of Lemma 2.2. We will derive (2.20). Relation (2.19) can be obtained in a similar way. Note that 

ESt(a) = trlog(l + F{a)) - logdet(l + e''), 



where 

Lemma 2.3 and simple algebra yield 



F{a) = e^/^e' 



k/2 a{kt-k) k/2 



Since 



we have 



^ESt(a) = tr (^(l + e-"('='-'^)e-'^') ^ (fc* - k) 



kt-k= / ^/ct„ du^t f e""''i[/i, fc]e~'*"'*dw, 







da 

Using ESf (0) = 0, integration yields Relation (2.20). □ 

2.3 Thermodynamic limit 

The proof of Proposition 1.3 is standard and we shall omit it (see [La]). To prove Proposition 1.4, we 
consider hnj as an operator on f) = £^(Z) {hi\j acts as zero on the orthogonal complement of t)A in f)). 
Clearly, Hm — > h strongly and so the strong limits 

s - lira ICp^t,M {a, u) /Cp,t(a,M), s - lim/CEs,t,M(a, ") = /CEs,t(a, 

AZ-i-oo AZ-i-oo 

exist. Moreover, /Cp 4(0;, u) and /CES,t(Q!, u) are given by the formulas (2.16), (2.17), (2.18) and hence are 
norm continuous functions of (p, t, a, u). Since i[knj, Hm] = i[k, h] does not depend on M and is a finite 
rank operator, an application of dominated convergence yields 

pa />! 

e#,t(a) = lim e^^t,M{a) ^ t d-f dutr {K.^^t{l,u)i[k,h]) , (2.21) 

where # stands for p or ES and eEs,t{ct) = ESt(a). The functions 6:^^4(0;) are jointly continuous in 
{p, t, a) and real-analytic in a. The Trotter product formula yields 

lim ep^t{a) = eoo,t(a), 

p— )-oo 
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and Parts (l)-(6) of Proposition 1.4 follow. To prove (7) and (8), note that for given # in {p, ES} and 
t there is an e > such that tjv/(a) extends analytically to the ball \a\ < e in C and satisfies 
supjv/^o |a|<e \^#-t,M < CO- This observation, Parts (6) and (7) of Proposition 1.2 and Vitali's conver- 
gence theorem (see, for example. Appendix B in [JOPP]) imply (7) and (8). To prove (9), we note that the 
function 

= det (1 + (1 + e-''^')-\e-°''""e-'*''^'e'^''^'e'*'"'' - 1)) , 

is entire analytic. The bound | dct(l + A)| < c'l'^lli, where denotes the trace norm of A, together 
with the formula 



(where v is finite rank) imply that for any bounded set B C C, 



sup 

aeB,M>0 



ge2,t,A/(a) 



< OO. 



By Vitali's convergence theorem the sequence of characteristic functions of the measures Pj^m converges 
locally uniformly towards an entire analytic function. The existence of the weak limit Ft follows (see 
e.g.. Corollary 1 to Theorem 26.3 in [Bi]) and the convergence of the moments is a direct consequence of 
Vitali's theorem. Finally, (10) is a general fact which follows from Araki's perturbation theory of modular 
structure, see [JOP] for the proof and additional information. We shall not make use of (10) in this paper. 



2.4 The Gallavotti-Cohen functional 

Note that 

GCt(a)= lim lim log m fe~"*^''>0 ■ 

Since 

det(l + gfe-s,A//2gQ;(fct,M — fcjvj)gfc-s,jvj/2^ 



Ws A/ 



det(l + e'=) 

the arguments of the last two sections yield 

^lim loga;,,,M (c""*^''*^) = -< d-f dutr ^(l + c~''('='(i-")-'=-*"'e-'''-<=+'")) ^i[k,h] 

The operator k is bounded, commutes with ho and Ran (fc) C f)ac(^o)- It follows that 

s - lim fc, = s - lim e'^'^e"'^'"' W''" e"'"'' = w± kw^ = k± , (2.22) 

s— >-±oo >dioo 

and the dominated convergence theorem yields 

GCt{a)=t dj dutr(/CGC,t(7,w)i[fc,/i]), (2.23) 
Jo Jo 

where 

/CGC,t(a,w) = - fl + e-"('='(i-")-'=-'")e-'^-) . (2.24) 
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2.5 Main results 

Proof of Theorem 1.6. It follows from (2.18), (2.24) and (2.22) that for u e]0, 1[, 

/C+(a) =s-lim/CEs.t(a,u) = s - lini/CGC,t(a, = - (l + e-"(*^+-'''-)e-''- 
Hence, by dominated convergence, 

e+(a) = lim -EStia) = lim -GCtia) = f tT{IC+ij)i[k,h])d-f. 

t— >00 t t-i-OO t Jq 

It follows from ( 1 . 1 0) and Theorem 4. 1 in [ AJPP2] that 

tr (/C+(7)i[fc, h]) = -tr (^(l + e-T(«*'=^-'=)e-'=) i[fc, h]w^^ 

= l^tr (^(l + e-^(-*(^^)^o(^^)«(^)-^-o(£;))e-'^-o(^))"' {s*iE)koiE)siE) - fco(i?))) 

Theorem 4. 1 in [AJPP2] is quite general and in the special case considered here the above formula can be 
also checked by an explicit computation [La]. By Lemma 2.3, 

tr ^(l + e-^("*(^)'=°(-^)''(^)-'^'''(^»e-'^°(^))~^ {s*iE)ko{E)siE) - koiE))j 

= iLtrlog + e'">^E)/2^jis'iE)ko{E}siE)-koiE))^koiE)/2'^ ^ 

and so, Fubini's theorem yields 

e+(a) = / tr(/C+(7)i[fc,/i])d7 
Jo 

^ trlog (l + c'">^E)/2^y{s'{E)ko{E)siE)-ko{E))^ko{E)/2'^ 



d7 

det (1 + e'=o(B)/2ga{s-(£;)feo(iJ)s(£;)-fco(B))efco(B)/2^ 



27r 



y^^"^ det(l + efeo(s)) 27r' 

The formula for ep.+ (a) is derived in the same way, starting with 

/Cp.+ (a) = s — lim ICp,t (a, w) 

i— f oo 

for p e]0, oo[, and 

t-i-oo V / 

This, together with the remark preceding Theorem 1.6, prove Part (1). 

The formulas derived in Part (1) clearly show that ep,+ and e+ are real analytic functions of a satisfying 
ep^+(0) = e+(0) = 0. As limits of convex functions they are also convex and hence have non-negative 
second derivatives. Thus, these derivatives either vanish identically or have an isolated set of real zeros. 
It follows that ep_+ and e+ are either linear or strictly convex . The symmetry ep,+ (a) — ep_+(l — 
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a), a consequence of Proposition 1.2 (4), implies ep,+ (l) = and hence ep,+ is strictly convex if it 
doesn't vanish identically. Regarding e^, an explicit calculation (see [La]) shows that if it doesn't vanish 
identically, then e" (a) > for all a e M. 

It follows from (2.21), (2.23) that there exists e > such that the functions e2,t(a)/t, ESt(a)/i and 
GCt{a)/t have analytic extensions to the complex disc |a| < e and are uniformly bounded on this disc for 
t > 0. Vitali's convergence theorem. Proposition 1.4 (7) and the fact that GCt(O) = — f(cr)_|- imply 

e;,+ (0) = lim ie;,,(0) = lim jl^S'M = eV(0) = lim jGC'M = + . 

Moreover, 

e^(0) = lim iES;'(0) = lim -GC';{Q), e'{+(0) = lim -e'lM, 

s-oo t i— soo t voo t 

and since 63 ^(0) = ES"(0) by Proposition 1.4 (8), we derive that 63^(0) = e" (0). Using that uj+ is 
T* -invariant, one easily derives 

-GC';{Q) = \ 1* ((a, - (a) + )(a - {a) + ))^ [l ds 

{{(Ju - (cr) + )(cr - (cr) + )) du 



ds. 



This completes the proof of Parts (2) and (3). 

To prove (4), we compare e+(l) and eoo,+ (l)- Since det(l + A) = 1 + tr(yl) + det(A) for any 2x2 
matrix A and det{Ka{E)) = det{Ka^oo{E)), we have, taking into account Remark 2 after Theorem 1.6, 

det (1 + K^{E)) - det (1 + K^^^{E)) = tr {^koiE) ^,.1:.{E)\^ _ |^g/=o(s)+aA(i5) 

with A(£;) = s*{E)ko{E)s{E) - fco(£'). By the Golden-Thompson inequality (see Corollary 2.3 and 
Exercise 2.8 in [JOPP]), the above difference of traces is strictly positive unless 

[ko{E),A{E)] = [ko{E),s*iE)ko{E)s{E)]^0. (2.25) 

One easily verifies that for E this happens if and only if s{E) is either diagonal or off-diagonal. It 
follows that 

e+(l) > eoo.+ (l) = 0, 

unless h is reflectionless. 
To prove Part (5), note that 

det(l + K^^p{E)) - det(l + K^^g{E)) = tr(i^„,p(£;)) - tr(A'„^,(^)). 

The Araki-Lieb-Thirring inequality ([Ar5, LT], see also Theorem 2.2. and Exercise 2.8 in [JOPP]) implies 
that the function 

]0, 00]3p^ t^K^AE)) = tr ^^il-c.)ko{E)/p^2^siE)ko{E)s'iE)/p^il-a)koiE)/p^ ^/"^ ^ 

is strictly decreasing unless (2.25) holds, in which case this function is constant. Part (5) follows. □ 

Proof of Corollary 1.7. Part (1) follows from Theorem 1.6 and the Gartner-EUis theorem (see Appendix 
A. 2 in [JOPP] or any book on large deviation theory). The analyticity argument used in the proof of Part 
(3) of Theorem 1.6 and Bryc lemma (see [Br] and Appendix A.4. in [JOPP]) yield Part (2). □ 
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